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ABSTRACT 


This report summarizes the work performed under NASA iLewis Grant 
Nomber NAG 3-67 entitled, "Prediction of Sound Radiation from Different 
Practical Jet Engine Inlets" during the two year period from June 2, 1980 to 
July 31, 1982. The objectives of this study included the following: (1) 
upgrade existing computer codes which were based upon a special integral 
representation of the external solutions of the Helmholtz equation, so that 
they would yield accurate results for the acoustic radiation patterns in the 
field surrounding an axisymmetric body for non-dimensional wave numbers, 
based on duct radius, of up to twenty; (2) check the accuracy of these 
computer programs by the use of the point source method for the generation 
of "exact" solutions and then by compeurison with the results of other 
experimental and theoretical studies; and, (3) use these computer codes in a 
parametric study of the dcin .'^dence of the radiated sound field on input 
modal distribution, wave number, and inlet lip shape. The results of this 
study show that: (1) as the wave number (i.e., cut-off ratio) is increased for 
a given input modal distribution that the acoustic radiation peak moves 
towards the inlet centerline and becomes more compact (i.e., narrows); (2) 
as the input mode number is increased for a given cut-off ratio the acoustic 
radiation peak moves away from the inlet centerline and becomes more 
compact; and, (3) as the inlet lip becomes thicker the acoustic radiation 
peak moves towards the centerline of the inlet and becomes more compact. 
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I. INTRODUCTION 


This report summarizes the work performed under NASA Lewis 
Contract Number NAG 3>67 during the period Tune 2, 1980 to July 31, 1982. 
This research program was concerned with the determination of the 
radiation patterns around different practical jet engine inlet configurations 
when subjected to different excitation conditions. This was accomplished in 
two phases, with the first phase completed during the first year of this 
contract and the second phase the following year. 

During the first year of this project the computer codes required for 
this investigation were developed. The initial task of this first phase was the 
upgrading of existing computer codes for the calculation of both the surface 
and field distributions of the acoustic quantities of interest (i.e., the 
acoustic potential and acoustic velocity normal to the wall). The theoretical 
basis for these programs is presented in Refs. 1-3. These references track 
the development of a unique integral formulation of the acoustic radiation 
problem and its application in the calculation of two dimensional, three 
dimensional and cylindrically symmetric radiation patterns. The cylindrically 
symmetric formulation of Ref. 3 (See Appendix A of Ref. tt.) was used in this 
investigation as all of the inlet geometries investigated were axisymmetric. 
The details of the upgrading of the computer codes is presented in Ref. 4. 
This initial computer code development work was accomplished using the 
NASA Lewis QCSEE inlet configuration of Ref. 5 ana was checked against 
"exact” solutions generated by the point source metiiod as described in Ref. 
3. 
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The second task completed during the first year was the consolidation 
of the two main computer programs, along with the many smaller computer 
programs developed to do such things as calculate the cut-off wave numbers 
for various modes, the generation of input modal distributions and the 
geometrical shapes of the bodies of Interest, and the data plotting routines 
into a coherent package of programs. This program package was then used to 
investigate the acoustic behavior of five different inlet configurations; 
namely: the NASA Langley Bellmouth^^^; the NASA Lewis 3T15D-1 Ground 
Test Nacelle ’ ; and three hyperbolic inlets of 50, 70 and 90 degree 
divergence angles, similar to those investigated by Y. C. Cho^®*^^. Although 
the main objective of this study was to compare the radiation patterns 
generated by the integral technique with thoi>e calculated/measured in 
other theoretical/experimental studies, reflection coefficients were also 
calculated at the duct exit plane for certcun specific inlet configurations 
where these results were available. The results of these studies are 
presented in Ref. 12. 

During the first half of the second contract year a series of 35 
computer runs were performed for the NASA Langley Bellmouth inlet 
configuration^’^. These 35 computer runs provide a parametric study of the 
dependence of the sound radiated from this inlet geometry upon the 
specified input modal distributions and the cut-off ratios given in Table I. 
For each of these runs, both the distribution of the reflection coefficient at 
the duct exit plane and the distribution of the SPL (Sound Pressure Level) in 
the field were calculated. These results were cross plotted holding the input 
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modal distribution or the cut-off ratio constant and can be found in Ref. 13. 

During the final half year of this two year contract the computer 
programs were updated to run on the new computer systems here at Georgia 
Tech as the old CDC Cyber 70/74 system was replaced by two new CDC 
Cyber computers; a 730 and a 760. Also, with the new campus network in 
operation, printouts could be obtained on site rather than only at the 
computer center which greatly reduced printing costs. This network also 
allowed the down loading of data files from the CDC computers to the 
Hewlett-Packard 2100 series mini-computer in the School of Aerospace 
Engineering for plotting. Having accomplished this, the computer program 
which uses the calculated surface distributions on the body to calculate the 
radiated sound field outside the body was modified at the sponsor's request 
to: (1) calculate the acoustic radiation pattern all around the body from 0 to 
180 degrees from the inlet centerline at one degree increments instead of 
from 0 to 90 degrees in two degree increments; and, (2) to calculate the 
radial distributions of the reflection coefficients across the inlet at five 
axial stations from the striaght duct exit plane to the inlet exit plane. This 
updated field code was then re-run for the 33 cases done previously for the 
NASA Langley Bellmouth inlet configuration. The field data was then 
plotted and cross plotted as before (See Ref. 13.) holding either the input 
modal distribution or the cut-off ratio constant. The distributions of the 
reflection coefficients at the five axial stations were also plotted for the 
five cases run for each of the two input modal distributions (i.e., (M, N) s 
(4,1) and (4,3)). 
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Two elliptical inlei geometries were supplied by the sponsor for study 
and they will be designated here as the "thin” and "thick" lipped elliptical 
inlets. Thirty five computer runs were done for each of these inlet 
configurations at the same input conditions as was done for the NASA 
Langley Bellmouth inlet configuration (See Table I.). The same plots and 
cross plots were obtained for each of these geometries and the results are 
presented herein. In addition^ the radiation patterns generated by all 103 
computer runs were cross plotted varying the inlet geometry while holding 
the input modal distributions and the cut-off ratios constant. These plots are 
also presented in this report. 


II. THE LANGLEY BELLMOUTH INLET 


Figures 1-3 represent the distributions of the calculational points on 
the surface of the theoretical model used to simulate the NASA Langley 

Bellmouth inlet for the various computer runs. The open circles are the 

( 

points where the acoustic quantities are calculated while the locations 
denoted by the small crosses are those used in the integration procedure. 
The arrows represent the outward normals from the body at the various 
points on the body as these plots offer a visual check of the geometric input 
data for the computer programs. One should note that the driver plane also 
represents the duct exit plane as there is no straight duct section connected 
to the inlet lip. This was done to conserve computer space, which was 
limited in the old Cyber computers, and computing time, which goes up 
roughly as the square of the number of points used on the surface of the 
body, as the effect of a straight duct connected to the inlet lip was not 
investigated in this study. Finally, it should be noted that the theoretical 
geometry is terminated in all cases by a 2:1 ellipse as the Helmholtz 
equation is elliptic and is therefore only applicable to closed bodies. 

Descriptions of the computer runs conducted for each of the three 
geometries investigated are presented in Table II. Since the required 
computer runs for this parametric study cover a wide range of input modes 
and cut-off ratios (i.e., wave numbers), the number of calculational points on 
the body was not the same for all of the runs. This is due to the fact that as 
the wave number (cut-off ratio) is increased more points must be used on die 
body in the longitudnal direction to adequately discribe the distributions of 


5 


the acoustic quantities on the ^rface of ttie body. Aiso» as the input mode 
number is increased, the number of f>oints in the tangential direction must 
be increased for the same reason. 

In Figs. 4 and 5 the modulus and phase of the radial distributions of 
the reflection coefficients are presented at five axial locations along the 
inlet starting at the duct exit plane, (i.e., the driver plane), station i, and 
ending at the inlet exit plane, station 3. These stations are separated by 
equal axial distances so that station 3 is located at the middle of the inlet 
lip. The non-dimensional axial distances of these stations from the driver 
plane are given in Table III where the reference distance is the radius at the 
driver plane. Figures 4 and 5 present the results for the input modes of (M,N) 
= ((^,1) and (M, N) = (4,3), respectively. These plots show that as the acoustic 
wave moves out of the inlet the modulus of the reflection coefficient 
steadily decreases as does the phase. Also, as the wave number (cut-off 
ratio) goes up the modulus of the reflection coefficient decreases while its 
phase increases. 

The radiation patterns for the 35 cases run for the Langley Bellmouth 
are presented in Ref. 13 (See Figs. 4-38.) and therefore will not be repeated 
here. When these data are displayed on separate plots the results are not 
easy to compare and therefore trends are not obvious. To alleviate this 
problem, these data are compared in plots of the relative SPL (Sound 
Pressure Level) in decibels and the phase in degrees versus the angular 
location in Figs. 6-10. Each of these plots describes SPL levels referenced to 
the point in the field where the maximum SPL occurs for each case. Each 
figure contains the SPL in dB plotted in two different ways so that trends 
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can easily be detected. For example, the same data are plotted in Figs. 6a dc 
c and 6b & d where in Figs. 6a & b the curve are all referenced to the same 
SPL level (i.e., zero) in Figs. 6c 6c. d each curve is referenced to a different 
SPL level. The actual data was calculated at one degree increments on a 
circle centered at the inlet exit plane with a radius of 20.33 driver radii. 
Each figure contains all of the results for the various cases run for each 
input mode (i.e., Fig. 6 contains the results for an input mode of (M,N) = 
(1,1), etc..). One should note that there is a change of scale in Figs. 6a & b 
and in Figs. 7a & b and that some of the cases are repeated for clarity as is 
done in Figs. 6c & d, 6e & f, 7c 6c d, and 7e Sc f. One should also note that 
since the theoretical geometry outside and behind the inlet lip does not 
model the true geometry of the actual inlet, the results in the field past 90 
degrees from the centerline are suspect and at best only provide the 
qualitative trends of the sound distribution in these locations. 

In Figs, ii-16 the relative SPL and phase distributions of the sound 
fields associated with the Langley Bellmouth are presented. Specifically, 
these figures describe cases where the cut-off ratio is held constant and the 
input modal distribution is varied. As in the previous set of figures, the 
relative SPL is plotted in two different ways in each figure for clarity. 
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III. THE THIN LIPPED ELLIPTICAL INLET 


O 


The model geometries used for the 35 cases run for the thin lipped 
elliptical inlet are given in Figs. 17-19. As with the Langley Bellmouth inlet, 
three different numerical approximations were employed for various input 
conditions. The runs for which each approximation was used are presented 
in Table II. One should note that the nominal point spacings on this inlet are 
half of those used on the Langley Bellmouth thereby making the calculations 
more accurate. This was done for the following reasons: (1) with the new 
computers more points could be taken on the surface of the body as more 
computational space was available; and, (2) since the inlet lip is 
considerably smaller, a closer point spacing was necessary to adequately 
describe the relatively faster changing curves. This inlet lip consists of a 2: 1 
ellipse whose semi-major axis is 0.2 of the non-dimensional driver radius, a. 
It is capped with a semi-circle of radius 0.025a and has a 2:1 ellipse as its 
outer wall and is terminated by yet another 2:1 ellipse. 

In Figs. 20-24 the SPL in dB are plotted and tabulated from 0 to 90 
degrees in two degree increments. Each figure contains all of the plots for 
the cases run for each modal input (e.g.. Fig. 20 consists of the 9 plots, one 
for each cut-off ratio, for an input mode of (M,N) ~ (1,1), etc..). These data 
represent actual values generated in each case and not just the relative 
values. Also, when comparing these plots, care must be exercized as all of 
the plots do not use the same limits. 

The modulus and phase of the radial distributions of the reflection 
coefficients are presented in Figs. 25 and 26 for input modal distributions of 
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(MfN) s (4,1) and (4,3), respectively. The results for this inlet configuration 
show that the reflection coefficients don't follow the same trends over the 
entire spectrum of investigated wave numbers for both input modal 
distributions. For example, the modulus of the reflection coefficient 
decreases as the wave moves down the inlet in Figs. 25a,c,e and g. while it 
increases in Figs. 25 i and k. 

In Figs. 27*31 the relative SPL in dB and phase in degrees are plotted. 
Each figure represents all of the cases run for the thin lipped elliptical inlet 
for a given input modai distribution. As before, the relative SPL is plotted in 
two ways In each figure for the sake of clarity. Also, as before, some of the 
curves are replotted for certain cases (e.g., see Fig. 27a & b) for ease of 
comparison and because the scales change between these plots. It is of 
interest to note that although the lobe which contains the maximum SPL 
changes for the Langley Bellmouth inlet (See Fig. 10a & b.) when the (M,N) = 
(4,4) mode is cut-on (See Table I.), this does not happen for the thin lipped 
elliptical inlet. While in the Langley Bellmouth inlet the maximum SPL 
moves from the first lobe (closest to the centerline) to the third lobe as the 
wave number increases, the maximum SPL always remains in the third lobe 
for the thin lipped elliptical inlet. 

The radiation patterns are re-plotted in Figs. 32-37 holding the cut- 
off ratio constant and varying the input mode. Again, the relative SPL is 
plotted in two different ways in each figure so that trends are more easily 
discernable. 
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IV. THE THICK LIPPED ELLIPTICAL INLET 

Figures contain the model of representations used for the thick 
lipped elliptical iniat geometry. As with the other two inlet geometries^ 
these three different numerical descriptions of the model of the geometry 
are used for certain of the 35 computer runs done for this geometry (See 
Table II.). The inlet lip is simply twice as large m the thin lipped inlet 
geometry so that It too is described by a 2:1 ellipse whose semi-major axis is 
O.^a instead of 0.2a. It is also capped with a semi-circle (of radius 0.05a) and 
terminated by a 2:1 ellipse. 

The SPL in dB is both poltted and tabulated for the 35 cases run for 
this geometry (See Table I.) in Figs. 41-<^5. As before, each figure contains 
the plots for all of the cases run for a particular input mode. One should 
note again that these plots are of the actual SPL values calculated in the 
programs and that the scales on all of the plots are not the same. 

In *?igs. and (^7 the modulus and phase of the radial distributions of 
the reflection coefficients at five stations down the inlet are presented. 
Each figure contains the plots of all of the cases run for a particular mode. 
As with the thin lipped elliptical Inlets results, these distributions do not 
follow consistent trends as the wave number is increased or as the waves 
move down the inlet. The axial locations of the different stations are given 
in Table III. 

The relative SPL in dB and phase in degrees of the acoustic radiation 
patterns outside the inlet are plotted in Figs. ^S-52. As before, each figure 
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represents all of the cases run for a given mode for the thick lipped 
elliptical inlet. Also, the relative SPL is plotted in two ways in each figure 
and in Pigs. 48 and 49 some of the curves are duplicated in separate plots 
for ease of comparison. 

In Figs. 53-58 the radiation patterns are re-plotted holding the cut- 
off ratio constant in each case rather than the input mode. Each figure 
consists of three plots; two of the relative SPL in dB and. one of the phase in 
d^rees. The data for these plots were generated at one degree increments 
and the plots were obtained by having the computer draw straight lines 
between these data points. This accounts for the sharp corners observed on 
some of the curves (See Fig. 58a.) and the gaps in the curves when they cross 
180 degrees (See Fig. 58c.). 
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V. THE EFFECT OF INLET LIP CURVATURE ON ACOUSTIC 
RADIATION PATTERNS 


The principle goal of this parametric study was to show the effect of 
inlet lip curvature on the acoustic radiation patterns generated by the same 
acoustic excitation.. In the course of this study 35 computer runs were made 
for each of the three inlet geometries; the Langley Bellmouth, the thin 
lipped elliptical inlet, and the thick lipped elliptical inlet; for the same input 
model distributions (See Table I.). Since the duct exit/driver plane was non- 
dimensionalized to one for each geometry the acoustic excitation applied in 
the computer program was the same for each geometry. This does not mean, 
however, that the same acoustic power was supplied in each case as this is 
an elliptic problem and conditions at the driver plane depend upon the 
properties on the remaining boundaries. In these studies the acoustic 
potential (which is directly related to the acoustic pressure) was 6pecified as 
the boundary condition at the driver plane while a zero normal acoustic 
velocity (i.e., a hard wall) was specified on the rest of the body. This means 
that in each case the normal acoustic velocity was calculated on the driver 
plane which, in turn, resulted in a different acoustic power out of the driver 
(and therefore the inlet) for each case. 

Two options for normalization of the radiation patterns were 
therefore available: (1) normalize with respect to the peak radiated acoustic 
pressure; or, (2) normalize the curves with respect to the total radiated 
power (which for a hard walled geometry is equivalent to the power out of 
the driver). Since this study was more interested in the shape of the 
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acoustic radiation patternsi than in the total power radiated, the first 
normalization scheme was used. The second normalization scheme would 
entail the calculation of the PWL (Sound Power Level) in the field rather 
than the SPL (Sound Pressure Level) and than an integration under the curve 
for the normalization factor. It had previously been found by the authors of 

It 

this report that in general the PWL and SPL curves are very similar and 
almost interchangeable in the field surrounding the inlet. This being the 
case, the effect of the second type of normalization on the curves can easily 
be determined by inspection; that is, by doing a rough estimate of the area 
under the curves. One must be careful however as the area of the sphere 
must be taken into account so that the middle must be more heavily 
weighted than the ends. 

The results of the 105 computer runs are presented in Figs. 59-63. 
Each figure contauns all of the plots for each input mode, one plot for each 
cut-off ratio (wave number). An examination of these figures reveals two 
trends; which are: (1) as the inlet lip gets larger the maximum radiation peak 
moves further towards the centerline (in front of the inlet); and (2) as the 
inlet lip thickens, the maximum radiation pead< becomes narrower and better 
defined. Both of these trends are consistent with the expectation that a 
larger inlet lip will better focus the radiated acoustic energy. 
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